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Abstract 


Analysis of temperature-sensitive (ts) mutant viruses is a classic method allowing 
researchers to identify genetic loci involved in viral replication and pathogenesis. Here, 
we report genetic analysis of a ts strain of mouse hepatitis virus (MHV), tsNC11, 
focusing on the role of mutations in the macrodomain and the papain-like protease 2 
(PLP2) domain of nonstructural protein 3, a component of the viral replication complex. 
Using MHV reverse genetics, we generated a series of mutant viruses to define the 
contribution of macrodomain- and PLP2-specific mutations to the ts phenotype. Viral 
replication kinetics and efficiency of plating analysis performed at permissive and non- 
permissive temperatures revealed that changes in the macrodomain alone were both 
necessary and sufficient for the ts phenotype. Interestingly, mutations in the PLP2 
domain were not responsible for the temperature sensitivity but did reduce the 
frequency of reversion of macrodomain mutants. Co-immunoprecipitation studies are 
consistent with an interaction between the macrodomain and PLP2. Expression studies 
of the macrodomain-PLP2 portion of nsp3 indicate that the ts mutations enhance the 
proteasome-mediated degradation of the protein. Furthermore, we found that during 
virus infection, the replicase proteins containing the MAC and PLP2 mutations were 
more rapidly degraded at the non-permissive temperature, as compared to the wild-type 
proteins. Importantly, we show that the macrodomain- and PLP2-mutant viruses trigger 
production of type | interferon in vitro and are attenuated in mice, further highlighting the 


importance of the macrodomain-PLP2 interplay in viral pathogenesis. 
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Importance 


Coronaviruses are emerging human and veterinary pathogens with pandemic potential. 
Despite the established and predicted threat these viruses pose to human health, there 
are currently no approved countermeasures to control these infections in humans. Viral 
macrodomains, enzymes that remove post-translational ADP-ribosylation of proteins, 
and viral multifunctional papain-like proteases, enzymes that cleave polyproteins and 
remove polyubiquitin chains via deubiquitinating (DUB) activity, are two important 
virulence factors. Here, we reveal an unanticipated interplay between the macrodomain 
and the PLP2 domain that is important for replication and antagonizing the host innate 
immune response. Targeting the interaction of these enzymes may provide new 


therapeutic opportunities to treat CoV disease. 


Introduction 


Coronaviruses (CoVs) are enveloped, positive-sense, single-stranded RNA 
viruses that primarily infect the respiratory or gastrointestinal tract. CoVs can emerge 
from an animal reservoir, such as bats, to infect a new species and cause epidemic or 
pandemic disease with high mortality. Recent emergence events exemplified by Severe 
Acute Respiratory Syndrome coronavirus (SARS-CoV) and Middle East Respiratory 
Syndrome coronavirus (MERS-CoV) in humans (1), and Swine Acute Diarrhea 
Syndrome coronavirus (SADS-CoV) in domestic pigs (2), have demonstrated how 
devastating these viruses can be within naive populations. To date, there are no 


approved antivirals or effective vaccines that protect humans from coronavirus 
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diseases. Therefore, identifying viral factors that contribute to pathogenesis and 
characterizing novel targets for therapeutic interventions are two important approaches 


to facilitate the development of effect vaccines and antivirals. 


The murine coronavirus, mouse hepatitis virus (MHV), is widely used as a model 
system to study coronavirus replication and pathogenesis in mice. The replication of the 
virus initiates with the engagement of the spike glycoprotein with a host cell receptor 
and the release of the positive-sense RNA into the cytoplasm of the cell. The large (~32 
Kb) viral genomic RNA is translated to produce two long polyproteins, pp1a and pp1ab, 
which are processed by viral proteases, including the papain-like proteases (PLP1 
and/or PLP2), and the 3C-like protease (3CLpro or Mpro), into 16 nonstructural proteins 
(nsp1-16, Figure 1A). To generate the viral replication complex, the coronaviral nsps 
sequester host endoplasmic reticulum (ER) to generate convoluted membranes and 
double-membrane vesicles (DMVs), which are the sites of viral RNA synthesis (3, 4). 
The viral replication complex generates a nested-set of dsRNA intermediates to 
produce copious amounts of mRNAs, which are then translated to produce the 
structural (spike, envelope, membrane and nucleocapsid) and virus-specific accessory 
proteins. The genomic RNA and structural proteins assemble in the ER-Golgi 
intermediate compartment to generate infectious virus particles that are released from 


the cell (5, 6). 


CoV replication induces profound rearrangement of the host ER, and generates 
viral dsRNA intermediates, processes that can be sensed by the host to activate the 
innate immune response. As a result, CoVs have evolved multiple strategies to 


counteract and delay activation of these host immune responses and establish an 
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environment amenable to virus replication. These strategies include: expressing 
species-specific accessory proteins as modulators of innate immune responses 
[reviewed in (7)]; encoding highly-conserved nonstructural proteins that serve as 
interferon antagonists (8-10) and sequestering viral RNA in DMVs (3, 4) to prevent 
detection by host pattern recognition receptors. A key component in the assembly of 
the DMVs is nsp3 (11, 12). To date, eleven distinct nsp3 domains have been identified 
using either bioinformatic approaches or enzymatic studies (Figure 1A). Here, we focus 
on two multifunctional components encoded within nsp3, the macrodomain and the 


papain-like protease 2. 


The region researchers now term the macrodomain was originally identified in 
the 1990s as a highly-conserved domain of unknown function, termed the X domain, 
contained within the replicase polyprotein of rubella virus, hepatitis E virus (HEV) and 
coronaviruses (13-15). Structural and biochemical studies revealed that the X domain 
exhibited structural similarity to the cellular histone MacroH2a and catalyzed 
measurable ADP-ribose-1"-phosphatase (ADRP) activity (16-18), although the 
functional significance of this enzymatic activity was unclear. Using reverse genetics to 
inactivate the catalytic site of the enzyme, researchers found that ADRP activity was not 
essential for CoV replication in cultured cells (19). However, further studies revealed 
that an ADRP-catalytic mutant virus was attenuated in mice (20), and that ADRP activity 
in SARS-CoV and human CoV-229E mediated resistance to antiviral interferon 
responses (21). These findings were consistent with an essential role for enzymatic 
activity in vivo; although, the target for the ADRP activity was still unclear. A 
breakthrough came in 2016 from a study revealing that the macrodomain of hepatitis E 
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virus acts as an ADP-ribose hydrolase (22). ADP-ribosylation is a known post- 
translational modification that regulates cellular activities (23); therefore, viral enzymes 
that reverse this process could interrupt host-cell signaling. For CoVs, nsp3 
macrodomain activity was shown to promote MHV-induced encephalitis (24) and 


increase virulence during SARS-CoV infection (25). 


Another highly-conserved enzyme contained within nsp3 is the papain-like 
protease 2 (PLP2). For MHV, PLP2 is responsible for processing the nsp3/4 junction 
using a highly-conserved LXGG/X cleavage site (26). Studies using SARS-CoV 
revealed that the single papain-like protease encoded on nsp3 (termed PLpro) cleaves 
all three sites at the amino-terminal end of the polyprotein (27). PLpro also functions as 
a deubiquitinating enzyme (DUB), capable of removing polyubiquitin chains from 
substrates (28, 29). Structural studies revealed that CoV PLpro/PLP2s are similar to 
cellular DUBs (30). Enzymatic analysis revealed that CoV PLpro/PLP2s are 
multifunctional with protease, deubiquitinating and delSGylating activity (30-33). The 
viral DUB activity has been implicated as a modulator of the innate immune response to 
viral infection (82, 34, 35), but the target(s) of the DUB activity have not yet been 
identified. Thus, both the PLP2 and macrodomains of nsp3 have been independently 


identified as contributors to coronavirus virulence. 


In this study, we characterized a temperature-sensitive MHV mutant virus 
containing mutations within both the macrodomain and PLP2 domain. We investigated 
the contribution of these mutations to the temperature-sensitive phenotype as well as 
the resulting effects on viral pathogenesis. The results presented here reveal a 


previously undescribed interplay between the macrodomain and PLP2 domain that 
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impacts replication, antagonizes the innate immune response, and contributes to viral 
pathogenesis. Modulating the macrodomain-PLP2 interaction may provide new 


opportunities for therapeutic intervention. 
Results 


Identifying mutations associated with a temperature-sensitive phenotype. Murine 
coronavirus strain tsNC11 was generated by chemical mutagenesis, plaque purified, 
and validated as a temperature-sensitive (ts) mutant defective in positive-sense RNA 
synthesis at non-permissive temperatures (36). Complementation analysis indicated 
that tsNC11 harbored mutation(s) in the ORF 1a region of the replicase polyprotein, but 
the specific mutations were unknown. To identify the nucleotide changes in tsNC11, we 
isolated the genomic RNA from viral supernatant, subjected it to deep sequencing then 
aligned the reads to the genomic sequence of MHV-A59 (GenBank accession 
#AY910861). In agreement with the complementation study by Schaad et al. (36), the 
sequence analysis revealed 7 non-synonymous substitutions in the ORF1a of tsNC11. 
These substitutions resulted in 7 amino acid changes: two in nsp2 (14V and T5431), four 
in nsp3, and one in nsp10 (P23S). The four mutations within nsp3 are distributed 
between the macrodomain (K532E and G554D) and the PLP2 domain (D1026N and 
D1071N) (Figure 1A). As noted above, previous studies documented the importance of 
the macrodomain and PLP2 domain in virus replication and disease; therefore, we 
focused our efforts on evaluating how these substitutions contributed to the ts 
phenotype, the stability of the phenotype, and the pathogenesis of the virus in mice. 

To evaluate the contributions of the macrodomain and PLP2 domain mutations to 


the ts phenotype, three mutant viruses were generated using the MHV-A59 reverse 
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genetics system (37). The first mutant virus, designated MACmut, contains the 
macrodomain mutations K532E and G554D. The second virus was engineered with the 
D1026N and D1071N mutations within the PLP2 domain and is designated PLP2mut. 
The third virus, MAC/PLP2mut, combines the mutations in the macro- and PLP2 
domains into one virus. In addition, an isogenic wild-type MHV (icWT) was used as a 
control. These viruses were recovered, plaque purified, and propagated in DBT cells at 
a permissive temperature of 32°C. Deep-sequencing results confirmed the incorporation 
of the desired nucleotide changes in nsp3 and revealed no additional amino acid 
changes within the ORF1 region. 

First, we evaluated the one-step growth curves of all 5 viruses (tsNC11, icWT 
and the 3 engineered mutants) at the permissive (32°C) and non-permissive (37°C and 
40°C) temperatures. As expected, icWT replicates to high titer at all three temperatures, 
whereas tsNC11 is impaired at both 37°C and 40°C, as reported by Schaad et al. (36) 
(Fig. 1B). Analysis of the three engineered mutants revealed that the two substitutions 
in the PLP2 domain were not sufficient to confer a temperature-sensitive phenotype, as 
the kinetics of replication mirrored those of the wild-type virus. In contrast, the MACmut 
virus exhibited reduced virus replication at 40°C, but was only slightly impaired at 37°C. 
The MAC/PLP2 mutant virus mirrored the kinetics of tsNC11 with impaired replication at 
both 37°C and 40°C, implicating the mutations in both the macro- and PLP2 domains as 


contributors to the temperature-sensitive phenotype of tsNC11. 


We also evaluated the plaque size and efficiency of plating (EOP) of the viruses 
at permissive and non-permissive temperatures. As expected, all viruses replicated to 
high titer and formed similarly-sized plaques at 32°C (Figure 2, upper panel). tsNC11 is 
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profoundly temperature-sensitive, with a low number of plaques detected at the 10" 
dilution plate incubated at the non-permissive temperature. The tsNC11 plaques that 
were detected at the non-permissive temperature exhibited a large-plaque phenotype 
suggesting that these viruses may be revertants. We found that the PLP2mut virus 
formed large plaques at 40°C, which is consistent with the results of the kinetic analysis 
and indicates that the mutations in the PLP2 domain are not sufficient to cause the ts 
phenotype. Analysis of the MACmut virus revealed a mixed population of small and 
large plaques at the non-permissive temperature, the majority of which displayed the 
small-plaque phenotype. The MAC/PLP2 mutant virus mirrored the plaque size and 
plating efficiency of tSNC11. We calculated the EOP values, which represent the ratio of 
viral titers obtained at 40?C and 32?C (Figure 2B). Again, both the PLP2mut and icWT 
viruses had similar titers at both temperatures, resulting in an EOP of ~1. In contrast, 
the MACmut virus exhibited titers that were significantly lower at 40°C compared to 
titers obtained at 32°C (EOP=10%). These results indicate that the MACmut virus, but 
not the PLP2mut or icWT viruses, has a defect in plaque formation at the non- 
permissive temperature. Taken together, these data demonstrate that the mutations in 
the macrodomain, but not those in the PLP2 domain, are the major determinants of the 
ts phenotype of tsNC11. Additionally, these results are consistent with a critical role of 
the macrodomain in viral replication. Interestingly, we found that the MAC/PLP2 mutant 
virus mirrored the plaque size and low reversion frequency of tsNC11 (EOP = 10°), 
supporting a role for the PLP2 domain as a genetic enhancer of the ts phenotype. A 
genetic enhancer, as defined by genetic studies of eukaryotic organisms, is a mutation 


in one gene that intensifies the phenotype caused by a mutation in another gene (38). 
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Evaluating revertants of the MACmut virus. While generating the MACmut virus, we 
noticed that, in addition to the majority population having the small-plaque phenotype, a 
subpopulation of large plaques were also present at 40°C. The large plaques 
consistently appeared even after several rounds of plaque purification of the small- 
plaque isolates. Therefore, we asked if the small plaques were formed by temperature- 
sensitive viruses, while the large plaques were due to revertant viruses. To address this 
question, we selectively isolated plaques with different sizes and propagated them at 
32°C to obtain viral stocks for subsequent analysis (Figure 3). We found that the small- 
plaque isolates recapitulated the phenotype of the parental MACmut virus: small 
plaques and similar EOP values (Figure 3A). In contrast, the large-plaque isolates 
exhibited a phenotype similar to icWT. Sequencing results of PCR amplicons, 
representing the region spanning the macro- and PLP2 domains, revealed that small- 
plaque isolates had no additional mutations in either the macrodomain or PLP2. In 
contrast, the large-plaque isolates had either a true reversion (D554-to-G), or harbored 
putative suppressive mutations located within the macrodomain, or the adjacent, 
downstream sequence (Figure 3B). Among seven large-plaque revertants, all 
maintained the K532E mutation, indicating it was not associated with the ts phenotype. 
Three isolates had the D554-to-G reversion, suggesting that it may be sufficient for the 
ts phenotype of the MACmut virus. We found that isolates 4-7 maintained the 
engineered mutations, but had also acquired additional, potentially suppressive, 
mutations in the downstream region. Together, these results indicate that altering the 
coding sequence of either the macrodomain or the downstream region is likely sufficient 


to revert or suppress the ts phenotype caused by the G554D mutation. We also 
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evaluated the MAC-PLP2 region of tsNC11 large-plaque revertant viruses and found 
that all three isolates had the D554-to-G reversion (Figure 3C), consistent with our 


findings with the MACmut revertants. 


Mutations in PLP2 enhance the ts phenotype by reducing reversion frequency. 
We determined that the macrodomain mutations are the major contributors to the ts 
phenotype; however, we noticed that the MACmut virus did not completely phenocopy 
tsNC11. We found that the replication of the MACmut virus was defective at 40°C, but 
not at 37°C. In addition, the MACmut virus exhibited a higher EOP value (10°) 
compared to that of tsNC11 (10%) (Fig 2B), indicating a relatively high reversion 
frequency. These data imply that mutations outside the macrodomain may enhance the 
ts phenotype by stabilizing the replication defect, thereby preventing reversion to the 
wild-type phenotype (38). Therefore, we asked if the addition of the PLP2 mutations 
observed in tsNC11 could enhance the ts phenotype of the MACmut virus and reduce 
reversion. We found that the MAC/PLP2mut virus exhibits a severe replication defect at 
both 37°C and 40°C (Figure 1B), and only replicated under permissive conditions, 
similar to tSNC11. The MAC/PLP2mut and tsNC11 viruses exhibited similar EOP values 
(710^) (Figure 2B). Of note, the low titer of the MAC/PLP2mut virus at 40°C indicates a 
low level of reversion to the wild-type phenotype, suggesting that the PLP2 mutations 
stabilize the MACmut virus. Taken together, these data demonstrate that while the 
PLP2 mutations are not sufficient to cause the ts phenotype, they act to enhance the ts 
phenotype caused by the mutation in the macrodomain. Enhancing phenotypes have 


been described for other coronavirus interacting proteins (39, 40), which motivated us to 
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determine if the enhancement phenotype we detected here is due to an interaction 


between the macro and PLP2 domains. 


Evaluating macrodomain interaction with the PLP2 domain. The structures of 
several domains of nsp3 have been solved individually [reviewed in (41)] or in 
combination (42). However, owing to the size and complexity this protein, the complete 
structure of nsp3 remains unsolved. The capacity of the PLP2 mutations to enhance the 
ts phenotype in the presence of the macrodomain mutations raises the possibility of 
domain-domain interaction between the macrodomain and PLP2. To test this 
hypothesis, we generated plasmids that express either an epitope-tagged macrodomain 
(HA-MAC) or PLP2 domain (PLP2-V5) (depicted in Figure 4A). When these plasmids 
were co-transfected into HEK-293T cells, the expression of both the macrodomain and 
PLP2 were detectable by the cognate epitope antibodies (Figure 4B). We detected HA- 
MAC in lysates immunoprecipitated with anti-V5, and inversely, PLP2-V5 was detected 
when HA-MAC was immunoprecipitated from the lysates. These results indicate that the 
ectopically expressed macrodomain associates with PLP2 in cell lysates, consistent 


with either a direct or indirect interaction. 


Mutations in the macrodomain and PLP2 domain affect protein stability. Because 
we found that mutation in the macrodomain (G554D) is the major ts determinant and the 
PLP2 mutations enhance the ts phenotype, we reasoned that these mutations might 
alter protein folding, thereby rendering the protein unstable and susceptible to 
proteasome-mediated degradation. To determine if the mutations in the macrodomain 
and/or the PLP2 domain alter protein stability, plasmid DNA expressing wild-type or 


mutant forms of MAC/PLP2 polypeptide (Figure 5A) were transfected into HEK-293T 
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cells. The cells were maintained at 37°C throughout the experiment. We added 
cycloheximide (CHX) at 16 hours post-transfection to block translation, and harvested 
cell lysates at the indicated times. The level of expressed proteins was determined by 
immunoblotting (Figure 5B and C). The MAC/PLP2 (WT) protein was maintained at 
levels comparable to those prior to treatment, up to 5 hours post-treatment with CHX. In 
contrast, we detected rapid reductions in the levels of all of the mutant forms of the 
protein. Addition of the proteasome inhibitor MG132 blocked the degradation of the 
proteins (Figures 5B and 5C). These results indicate that mutations in both the MAC 
and PLP2 domains affect the protein folding and stability, rendering the proteins more 


susceptible to proteasome-mediated degradation. 


To determine if these MAC/PLP2 mutations affect the stability of the replicase 
proteins during virus replication at the non-permissive temperature, we performed 
temperature shift experiments as outlined in Figure 6. We infected cells with either WT 
or MAC/PLP2mut virus and incubated at the permissive temperature for 9.5 h. At this 
point, we added CHX to block translation and shifted the infected-cells to the non- 
permissive temperature. Cell lysates were collected every 30 minutes and evaluated 
using immunoblotting for the level of nonstructural intermediate nsp2-3 and product 
nsp3. We found that WT nsp2-3 and nsp3 were relatively stable, with loss of detection 
occurring at 3 hours after the temperature shift and addition of CHX (Fig 6C, lanes 2-8). 
In contrast, the levels of nsp2-3 and nsp3 in the MAC/PLP2mut-infected cells 
diminished more rapidly, with reduced levels at 1.5 hours after the temperature shift and 


addition of CHX (Fig 6C, lanes 9-15). These results support the finding that the MAC 
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and PLP2 mutations destabilize the replicase protein at the non-permissive 


temperature. 


ts mutant viruses induce interferon in macrophages and are attenuated in mice. 
Previous studies have shown that the papain-like protease domains of MHV, SARS- 
CoV, and MERS-CoV antagonize the IFN response, likely through the deubiquitinating 
activity of these enzymes (32, 34, 35, 43-46). In addition, coronaviral macrodomains 
have been shown to suppress IFN production both in vitro and in vivo (20, 21, 24, 25). 
We asked if the mutations in the macrodomain and PLP2 modulate the type | IFN 
response during infection of macrophages. As shown in Figure 7A, infection of mouse 
bone marrow-derived macrophages (BMDMs) with mutant viruses at permissive 
temperature produced significantly more IFN-a during infection compared to the icWT 
virus infection. At 12 hours post-infection, the MACmut virus induced 2-fold more IFN-a 
than icWT virus. Furthermore, the level of N gene transcript, which reveals the 
abundances of all viral mRNAs, was reduced in the MACmut-infected cells, compared 
to the wild type virus. We found that the PLP2mut virus elicited dramatically more IFN-a 
than WT virus, while the level of N gene expression was similar. The MAC/PLP2mut 
virus exhibited the most robust IFN-o and the lowest level of N gene expression. These 
results show that mutations in the macrodomain and the PLP2 domain result in elevated 
levels of type | IFN mRNA expression during infection of macrophages, further 


supporting the role of macrodomain and PLP2 in modulating host innate immunity. 


Because the ts mutant viruses had reduced replication efficiency and elicited 
type | IFN production during infection of macrophages, we were interested in evaluating 


the pathogenicity of these viruses. To this end, C57BL/6 mice were intracranially 
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inoculated with 600 plaque-forming units (PFUs) of virus and monitored for weight loss 
and mortality. As shown in Figure 7B, all WT virus-infected mice lost weight rapidly and 
succumbed to infection by day 11 post-infection. In contrast, the mutant virus-infected 
mice exhibited transient or no weight loss during the infection period and all mice 
survived. These results demonstrate that the ts mutant viruses are attenuated in vivo, 
and those mutations adjacent to the catalytic sites of the macro- and PLP2 domains can 


modulate viral pathogenesis. 


Discussion 


Identifying viral factors that modulate the immune response to viral infection 
provides new opportunities for developing novel antiviral interventions. Here, we 
described an unanticipated interplay between two previously characterized virulence 
factors, the macrodomain and the papain-like protease, of coronaviruses. The 
enzymatic activities of these domains have been implicated in removing post- 
translational modifications: macrodomains remove mono- or poly-ADP-ribose from 
proteins (18, 22, 23); deubiquitinating activity of viral papain-like proteases removes 
mono- or poly-ubiquitin chains from signaling proteins (30, 32, 47). Our study stems 
from characterizing a temperature-sensitive mutant virus that harbored mutations within 
both the macrodomain and the PLP2 domain of nsp3. We found that the mutation within 
the macrodomain (G554D) was associated with the most significant temperature- 
sensitive phenotype, but that this alteration of the macrodomain reverted to the wild- 
type phenotype at high frequency. However, viruses containing mutations in both the 


macrodomain and PLP2 domain reverted less frequently, consistent with the PLP2 
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domain having an enhancing effect on the ts phenotype. Although these two enzymes 
reside within the same nsp3 polypeptide (Figure 1A), to our knowledge, this is the first 
suggestion of an interplay between these domains. By expressing the macrodomain 
and papain-like protease 2 domain on independent expression plasmids, we were able 
to evaluate and detect co-immunoprecipitation of the proteins, consistent with either a 
direct or indirect interaction. Furthermore, we report that the mutations identified in the 
macrodomain and PLP2 domain destabilize the proteins, as revealed by proteasome- 
dependent degradation. Lastly, we demonstrate that these mutant viruses promote type 
| IFN production from macrophages in tissue culture and are attenuated in mice. This 
work confirms and extends previous studies that independently identified the 
macrodomain and the papain-like protease 2 domain as modulators of the innate 


immune response and virulence factors [reviewed in (31, 48, 49), (35)]. 


Macrodomains have been shown to play a role in the virulence of positive-sense 
RNA viruses including hepatitis E virus (HEV), alphaviruses, and coronaviruses 
[reviewed in (48, 49)]. Studies of the alphavirus Chikungunya virus (CHIKV) revealed 
that the macrodomain at the N terminus of nsP3 hydrolyzes ADP-ribose groups from 
mono-ribosylated proteins and that this de-ribosylating activity is critical for CHIKV 
replication in vertebrate and insect cells, and for virulence in mice (50). Interestingly, 
viruses engineered to encode a mutation of the CHIKV macrodomain catalytic site 
rapidly reverted to the wild-type sequence (51), similar to the high frequency reversion 
we reported for the MHV MACmut virus (Figure 3). Studies of the role of the 
macrodomain during coronavirus replication indicate that catalytic activity is not required 
for virus replication in interferon non-responsive cell lines (19, 20, 24). However, 
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catalytic activity is important for replication in primary cells and in mice, implicating the 
macrodomain in evading the innate immune response and promoting viral pathogenesis 
(20, 24, 25). Identifying the ribosylated substrates that are targeted by the viral 


enzymatic activity is an important future direction for this work. 


Our study implicated an adjacent viral domain, the papain-like protein 2 domain, 
as an interacting partner with the macrodomain. Interestingly, the helicase domain 
adjacent to the macrodomain of hepatitis E virus (HEV) was found to modulate 
macrodomain activity. Biochemical assays revealed that the presence of the HEV 
helicase domain in cis enhanced the binding of the macrodomain to ADP-ribose and 
stimulated the hydrolase activity (22). Furthermore, we previously found that the 
mutations in the Ubl-2 domain could cause a ts phenotype and destabilize the PLP2 
domain (52). Here, we found that the mutations in the macro- and PLP2 domains 
destabilized the replicase proteins, as shown by the more rapid degradation of the 
proteins after temperature shift. We speculate that there may be a dynamic interaction 


between adjacent domains within the nsp3 polyproteins. 


As a multidomain protein, nsp3 must hold a sophisticated architecture to function 
properly and precisely. To date, four essential functions have been documented for this 
multidomain protein: 1) interaction of the Ubl-1 domain with the nucleocapsid (N) protein 
is important for genomic RNA synthesis and encapsidation (39, 40); 2) proteolytic 
processing of the N-terminal region of pp1a and pp1ab to release nsp1, nsp2, and nsp3 
(26, 27); 3) hijacking the cellular reticular network in concert with other membrane- 
associated proteins (nsp4 and nsp6) to form virus-specific membrane structures for 


RNA synthesis (11, 12); and 4) antagonizing the innate immune response through the 
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405 


actions of the de-ADP-ribosylating activity of the macrodomain and the deubiquitinating 
activity of the PLP2 domain [reviewed in (31, 48)]. The removal of post-translational 
modifications such as ADP-ribosylation and poly-ubiquitination could be directed either 
at cellular proteins to redirect them for use during viral replication, or to subvert 
signaling of innate immune responses. Ultimately, structural and biochemical studies 
will be needed to fully investigate the multiple cis and trans interactions of nsp3 and to 
determine if there is a dynamic interplay that modulates the stability, substrate 


specificity and/or affinity of the enzymes and substrates. 


We found that the MAC/PLP2mut virus recapitulated the ts phenotype of tsNC11 
(Figure 1B). However, it is possible that some or all of the other mutations we identified 
by deep sequencing (14V and T5431 in nsp2 and P23S in nsp10) may contribute in a 
subtle way to the phenotype of tsNC11. Nsp2 was shown to be dispensable for MHV 
and SARS-CoV replication, but the deletion of the nsp2 coding sequence resulted in 
decreased viral replication and RNA synthesis (53). For nsp10, previous studies 
revealed that this protein plays critical roles in the 3C-like protease-mediated 
polyprotein processing and viral RNA synthesis (54, 55). The results from these studies 
indicate that the mutations in nsp2 and nsp10 may also contribute to a ts phenotype. 
While our study focuses on the contribution of the macrodomain and PLP2, further 
studies are needed to fully evaluate the impact of other ORF1a mutations on the 


replication and pathogenesis of coronaviruses. 


In summary, we report what is, to our knowledge, the first indication of an 
interplay between the macrodomain and papain-like protease 2 domain of CoV nsp3. 


We found that this interplay impacts virus replication efficiency, innate immune 
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antagonism and virulence in mice. A detailed understanding of the relationship between 
the macro- and PLP2 domains will require further structural and enzymatic studies. We 
anticipate that the genetic analysis, co-immunoprecipitation and in vivo pathogenesis 


outcomes reported here will facilitate these future studies. 


Materials and Methods 


Virus and cells. Human embryonic kidney (HEK) 293T cells (CRL-11268, ATCC) were 
cultured in Dulbecco’s modified Eagle medium (DMEM) with 10% fetal bovine serum 
(FBS) and 2% L-glutamine. Delayed brain tumor (DBT) cells were grown in minimal 
essential media (MEM) (catalog no. 21800-0400; Gibco) supplemented with 10% 
tryptose phosphate broth (TPB) media, 5% heat-inactivated FBS (Atlanta Biological), 
2% penicillin/streptomycin (pen/strep; HyClone), and 2% L-glutamine. The BHK-MHVR 
cell line was kindly provided by Mark Denison at Vanderbilt University Medical Center 
and cultured in Dulbecco’s modified Eagle medium (DMEM) (catalog no. 12100-046, 
Gibco) supplemented with 10% heat-inactivated FBS and G418 (0.8 mg/mL; HyClone). 
Differentiated BMDMs were maintained in bone marrow macrophage media containing 
DMEM (catalog no. 10-017-CV, Corning) supplemented with 30% L929 cell 
supernatant, 20% FBS, 1% L-glutamine, 1% sodium pyruvate, and 1% pen/strep. HeLa- 
MHVR cells (56) were grown in DMEM (catalog no. 12100-046, Gibco) supplemented 
with 10% FBS, 1% L-glutamine, 0.5% HEPES, and 1% pen/strep. Temperature- 
sensitive MHV strain tsNC11 was propagated in DBT cells at 32°C. The infectious clone 
MHV-A59 strain (GenBank accession no. AY910861) serves as wild-type (icWT) virus 


for this study. 
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Deep sequencing and bioinformatic analysis. Viral RNA was extracted from the 
supernatant of tsNC11-infected DBT cells incubated at 32°C. Isolated RNA was sent to 
GENEWIZ, Inc. for cDNA library preparation and Illumina Miseq high-throughput 
sequencing. Raw reads were subject to pairing and trimming and aligned to the genome 
sequence of the synthetic construct of MHV A59 strain (GenBank accession no. 
AY910861) using Geneious software (Geneious R7, https://www.geneious.com). A 
medium-low sensitivity and an iteration of up to 5 times were chosen. A total of 195,824 
sequences with a mean coverage of 898.8 were aligned to the MHV Synthetic Construct 
template. Polymorphisms were detected using the "find variations/SNPs" tool. 
Parameters included a minimum coverage of 5 with a minimum variant frequency of 
25% in order for a variation to be called. The maximum variant p-value was set at 10° 
and the minimum strand bias p-value was set at 10° when exceeding 65% bias. We 
focused our analysis on the first 12 kb of the replicase gene, since previous studies 


reported the ts phenotype was associated with changes in this region (36, 57). 


Generation of mutant viruses. All infectious clones were generated using the reverse 
genetics system previously established for MHV-A59 (37). Mutations identified by deep 
sequencing within the macrodomain and PLP2 domain were introduced into plasmids A 
and B, respectively, then verified by sequencing of the plasmid DNA. DNA fragments 
were ligated together and used for in vitro transcription of viral RNA. /n vitro transcribed 
genomic RNA and N gene RNA was electroporated into BHK-MHVR cells, which were 
overlaid onto DBT cells in a T75 flask. These cells were incubated at the permissive 
temperature of 32°C to facilitate the replication of ts mutant viruses. Supernatants were 
collected at the time when cytopathic effect was evident, usually between 36 and 48 h 
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post-electroporation. All infectious clone mutant viruses were plaque purified, 
propagated on DBT cells, and subjected to full genome sequencing to validate the 
genotype. These infectious clones were designated MACmut, PLP2mut, and 


MAC/PLP2, according the locations of introduced mutations as shown in Figure 1. 


Temperature-sensitive assay and one-step growth kinetics. To determine the 
temperature sensitivity of mutant viruses, the efficiency of plating (EOP = titer 40?C/ titer 
32°C) of virus was measured. DBT cells were seeded into two 6-well plates at 5.0 x 10° 
cells/well a day prior to infection. Each viral stock supernatant was serially diluted and 
inoculated onto the DBT cells. After 1 h incubation at 37°C, inoculum was removed, and 
cells were subsequently overlaid with 0.8% 2x MEM/agar mixture. One plate was 
incubated at 32°C for 60 h, and the second plate was incubated at 40°C for 48 h. 
Agarose-covered cells were fixed using 4% formaldehyde for 1 h and stained using 
0.196 crystal violet solution after removal of agarose. Plaques were counted and titers 


were calculated. 


To evaluate the kinetics of virus replication, a one-step growth curve was generated at 
each temperature. Briefly, DBT cells were infected with the designated virus at 
multiplicity of infection of 5 for 1 h at 37?C, then plates were incubated at the specified 
temperatures. The supernatants were collected at indicated time points and titrated on 


DBT cells incubated at 32?C for 60 h. 


Isolation and characterizations of ts revertants. To isolate ts revertants, plaque 
assays were performed at 40?C. Viruses from single plaques were isolated and 
propagated in DBT cells at 32?C to obtain viral stocks. To determine the ts phenotype of 


the isolates, the isolates were titrated at both 32°C and 40°C and the EOP values were 
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calculated as described above. To identify any mutations within the region of macro- 
and PLP2 domains of the revertants, viral genomic RNA was extracted using 
TriReagent (MRC, Inc.) according to the manufacturer's instruction and subsequently 
subject to cDNA synthesis. A genomic region (3976-6101 nt) containing the macro- and 
PLP2 domain was amplified by PCR using specific primers (Sense: 5- CAA GAA AGG 
TCT TTA GGG CTG CTT -3’; anti-sense: 5'- GAC ACC ATC AAC CTT CTC AAA TG - 
3’). The PCR products were sequenced and the sequencing results were compared to 


the tsNC11 sequence. 


MAC and PLP2 expression plasmids. Nucleotide sequences encoding the 
macrodomain [467-622 amino acids (aa) of nsp3] were amplified from a codon- 
optimized MHV nsp3 gene (sequence available upon request) and cloned into pCAGGS 
vector with an HA epitope tag, designated as HA-MAC. The pCAGGS-PLP2 plasmid 
(PLP2-V5) was generated in a previous study (52). The coding sequence of the 
macrodomain through PLP2 domain (467-1085 aa) was inserted into pcDNA3.1 and 
fused with a c-terminal V5 epitope tag (pMP-WT). Mutations were introduced into these 
constructs using site-directed mutagenesis PCR or Gibson Assembly technique to 
generate pMP-GD (G554D), pMP-2DN (D1026N/D1071N) and pMP-GD/2DN 


(G554D/D1026N/D1071N), which all contain a c-terminal V5 tag. 


Co-immunoprecipitation. HA-MAC and PLP2-V5 plasmids were co-transfected into 
HEK-293T cells in 35 mm dishes. Cells were harvested using 500 uL lysis buffer (20 
mM Tris pH 7.5, 150 mM NaCl, 1 mM EGTA, 1 mM EDTA, 1% Triton X-100, 2.5 mM 
sodium pyrophosphate, 1 mM f-glycerophosphate, 1 mM sodium ortho-vanadate, 1 


yug/mL leupeptin, 1 mM PSMF) and 200 ug of whole cell lysates were used for 
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immunoprecipitation with 1 ug primary mouse anti-V5 (R96025, Invitrogen) or -HA 
(MMS-101R-200, Biolegend) monoclonal antibody (Ab). Protein-Ab mixtures were 
rotated at 4°C overnight and then added 15 uL magnetic protein G beads 
(LSKMAGAO2, Millipore) for 1 h incubation. Beads were washed three times with 
washing buffer (lysis buffer composition except 450 mM NaCl) and eluted with 40 uL 2x 
sample buffer (10% glycerol, 5% B-ME, 3% SDS, 7.5 mg/mL Trizma-base, bromophenol 
blue). Eluted products and 5% of cell lysates as input were subject to SDS-PAGE gel 


electrophoresis and immune-blotting with anti-V5 or anti-HA antibodies. 


Evaluating protein stability after addition of cycloheximide. To determine the 
steady-state level of protein, 0.5 ug of the specified plasmid DNA was transfected into 
HEK-293T cells with transfection reagent TransIT-LT1 (MIR2300, Mirus) according to 
the manufacture’s recommendation. At 16 h post-transfection, cells were treated with 20 
ug/mL of cycloheximide (CHX) (5087390001, Sigma Aldrich) or a combination of 20 
yug/mL CHX and 10 uM MG132 (474790, Calbiochem), a proteasome inhibitor, and 
harvested at the indicated time points. Equal amounts of cell lysate were subjected to 
immunoblotting with anti-V5 or anti-B-actin (A00702, Genscript) antibodies. The relative 
intensity of MAC/PLP2 bands (relative to B-actin) were measured and calculated with 
AlphaView software (Protein Simple). To assess the rate of decay of the protein amount 
over time for the four viral protein types, we fit the two-parameter simple exponential 
nonlinear regression function, y = 0,e7%*, using the NLIN procedure in SAS 9.4 
software package and verified using Minitab software version 18. In this regression 


equation, 6, is the initial viral amount parameter at time zero and 6; is the slope or rate 
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of decay parameter. The slope parameters (02) were each compared with the WT slope 


using NLIN’s provided two-sided t-tests and p-values < 0.05 were deemed significant. 


To evaluate the steady level of replicase proteins nsp2-3 and nsp3, we 
performed a temperature shift experiment. Briefly, HeLa-MHVR cells were infected with 
either WT or MAC/PLP2mut virus (moi = 5) and incubated at the permissive 
temperature for 9.5 h, when the cells were shifted to 40°C and treated with 20 ug/mL of 
CHX. Whole cell lysates were prepared at 30 min intervals by the addition of lysis 
buffer A (4% SDS, 3% DTT, 40 % glycerol and 0.065 M Tris, pH 6.8). The lysates were 
passed through a 25-gage needle to break up aggregates, incubated at 37°C for 30 min 
and loaded onto a 6% SDS-PAGE, followed by transfer to a nylon membrane. The 
membrane was incubated with a 1:2,000 dilution of rabbit polyclonal anti-nsp2-3 
antibody (anti-D3) (58), followed by horseradish peroxidase (HRP)-conjugated donkey 
anti-rabbit IgG (H+L) (Southernbiotech), and developed with Western Lightening Plus- 
ECL reagents (PerkinElmer). The membrane was stripped and re-probed using a 
1:2,000 dilution of mouse anti-calnexin antibody, followed by HRP-conjugated goat anti- 


mouse IgG (H+L), and developed as above. 


Reverse transcription quantitative PCR (RT-qPCR). The protocol of RT-qPCR was 
described previously (9) with slight modification. Briefly, BMDMs were mock-infected or 
infected with wild-type or mutant MHVs at a multiplicity of infection (MOI) of 1 and 
incubated at a permissive temperature of 32 ?C. At indicated time points, cells were 
harvested for RNA extraction using an RNeasy Mini Kit (74104, Qiagen). An equal 
amount of RNA was used for cDNA synthesis using Rt2 HT First Strand Kit (330401, 


Qiagen). To determine IFN-a11, B-actin, or MHV-A59 N gene mRNA production, qPCR 
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was performed with specific primers for mouse IFN-a11 (PPMO3050B-200, Qiagen), 
mouse B-actin (PPM02945B-200, Qiagen) or MHV-A59 N gene (Sense: 5- AGC AGA 
CTG CAA CTA CTC AAC CCA ACT C -3'; anti-sense: 5- GCA ATA GGC ACT CCT 
TGT CCT TCT GCA -3” using RT2 SYBR Green qPCR Mastermix (330502, Qiagen) in 
the Bio-Rad CFX96 system. The thermocycler was set as follows: one step at 95 °C (10 
min), 40 cycles of 95 *C (15 s), 60 *C (1 min) and plate read, one step at 95 *C (10 s), 
and a melt curve from 65 °C to 95 °C at increments of 0.5 °C/0.05 s. Samples were 
evaluated in triplicate and data are representative of three independent experiments. 
The levels of mRNA were relative to B-actin mRNA and expressed as 24°" [ACT = 


Cr(gene of interest) T Cr(g.actin)]. 


Evaluating viral pathogenesis. The protocol for evaluating pathogenesis of MHV was 
approved by the Loyola University Chicago IACUC and previously described (59). 
Briefly, six-week-old C57BL/6 female mice were purchased from the Jackson 
Laboratory. Mice were intracranially inoculated with 600 PFU in 20 uL PBS and 
monitored daily for changes in body weight. Infected mice were euthanized when weight 
loss was over 25% according to the protocol. Statistical analysis of survival rate was 


evaluated using the log-rank test. 
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Figure legends 


Figure 1. Evaluating the replication kinetics of coronavirus temperature-sensitive 


mutants at permissive and non-permissive temperatures. (A) Schematic diagram of 


the MHV genome and the domains of nsp3. Abbreviations: Ubl1, ubiquitin-like domain 1; 


Ac, acidic region; PLP1, papain-like protease 1; MAC, Macrodomain; DPUP, domain 
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proceeding Ubl2 and PLP2; Ubl2, ubiquitin-like domain 2; PLP2, papain-like protease 2; 
NAB, nucleic acid-binding domain; G2M, coronavirus group 2 marker domain; TMDs, 
transmembrane domains; Y, coronavirus highly-conserved domain. Representative 
structures of the macrodomain with ribose (229E; PDB: 3EWR) and PLP2 (MHV; PDB: 
4YPT) are shown in cyan and green with catalytic pockets circled and the residues 
involved in catalysis shown in magenta. The mutations described in this study are 
shown in red. (B) Growth kinetics of MHV and mutants at three temperatures. DBT cells 
were inoculated with the indicated virus (MOI of 5) for 1 h at 37°C and then shifted to 
the indicated temperatures. Culture supernatant was collected at the indicated hours 
post-infection and titrated in DBT cells at 32°C. The data are representative of two 


independent experiments. Error bars indicate +SD. 


Figure 2. Analysis of plaque size and efficiency of plating at the permissive and 
non-permissive temperatures. (A) Representative plaque assays at 32°C and 40°C 
for icWT, tsNC11 and engineered mutant viruses. The dilution of the viral stock is 
indicated and selected to visualize ~20-50 plaques per plate. (B) Efficiency of plating 


(EOP) = average titer at 40°C/ average titer at 32°C. 


Figure 3. Analysis of small- and large-plaque variants in the MACmut virus 
population. (A) MACmut isolates with distinct plaque sizes were evaluated for a ts 


phenotype. Sequence analysis of individual plaque-purified revertant isolates identified 
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mutations in the macrodomain and the adjacent downstream region in the large-plaque 


variants of the MACmut (B) and tsNC11(C) viruses. 


Figure 4. Evaluating co-immunoprecipitation of the macrodomain and the PLP2 
domain. (A) Schematic diagram of the individual constructs used to evaluate potential 
interactions between the macrodomain and PLP2. (B) Western blotting to identify 
expression and co-immunoprecipitation of HA-MAC and PLP2-V5. HEK-293T cells were 
transfected with the indicated plasmid DNAs, lysates were prepared at 18 hours post- 
transfection, subjected to immunoprecipitation with the indicated antibody and the 
products analyzed by SDS-PAGE and immunoblotting. The data represent the results of 
three independent experiments. Astersks indicate the cross detection of IgG chains by 


secondary antibody. 


Figure 5. Mutations in macrodomain and PLP2 enhance degradation of the 
polypeptide. (A) Schematic diagram of constructs used to evaluate protein stability. (B 
and C) Western blotting detecting wild-type or mutant forms of MAC-PLP2 polypeptide 
in the presence of cycloheximide (CHX) or a combination of CHX and a proteasome 
inhibitor MG132. HEK-293T cells were transfected with the indicated expression 
plasmid of wild-type (WT) or mutant forms of MAC-PLP2. At 16 h post-transfection, cells 
were treated with 20 ug/mL of CHX or a combination of 20 ug/mL CHX and 10 uM 
MG132 and harvested at the indicated time points. Equal amount of cell lysate were 


subjected to immunoblotting with anti-V5 or anti-B-actin antibodies. The relative intensity 
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of MAC/PLP2 bands (relative to B-actin) were measured and calculated with AlphaView 
software. The experiment was repeated two times and the representative immunoblots 
(B) and the curves of relative intensity (C) are shown. The slope parameters of the 
decay curves were evaluated using non-linear regression and two-sided t-tests 


compared to WT. **, P<0.005; ****, P<0.0001. 


Figure 6. Mutations in the macrodomain and PLP2 alter the stability of replicase 
protein nsp3. HeLa-MHVR cells were infected with either icMHV-WT or 
MAC/PLP2mut virus (MOI of 5) and incubated at 32°C for 9.5 h, then 20 ug/mL of 
cycloheximide (CHX) was added and cells were shifted to the non-permissive 
temperature. Lysates were prepared every 30 min, and the proteins separated by SDS- 
PAGE, and nonstructural proteins nsp2-3 and nsp3 were visualized by immunoblotting. 
A) Schematic diagram of MHV replicase polyprotein indicating the processing pathway 
and the region identified by the anti-nsp2-3 antibody. B) Outline of the experiment. C) 
Western blot evaluating the level of nsp2-3 and nsp3 proteins detected after shift to the 
non-permissive temperature. This is representative data of two independent 
experiments. Arrowhead indicates detections of cellular protein in all lysates. Asterisk 
indicates degradation products detected by anti-nsp2-3 antibody in the MAC/PLP2mut 


virus-infected cells. 


Figure 7. Macrodomain mutant viruses induce type | interferon in primary 


macrophages and are attenuated in mice. (A) Mouse bone marrow-derived 
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macrophages were infected with the indicated virus (MOI of 1) at 32°C. Total RNA was 
extracted at the indicated time points and subjected to RT-qPCR. The mRNA levels of 
IFN-a (left) and N gene (right) are presented relative to B-actin. The results are 
representative of three independent experiments and subjected to a two-tailed, unpaired 
t-test. Error bars indicate +SD. ***, P<0.001; ****, P<0.0001. n.s.: not significant. N.D.: 
not detected. (B) Six-week-old mice were injected intracranially with either icWT or the 
indicated ts mutant virus (600 PFU per mouse) and monitored for weight loss. Viral 
pathogenicity was evaluated by body weight loss (left) and percent survival (right). The 
number (n) of infected mice is indicated in parentheses. Error bars indicate + SEM. 


Differences in survival rates were calculated using a log-rank test. 
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